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1 6  SUPPLEMENTARY  NOTATION 


I  1 7  COSATI  CODES 

FIELD 

GROUP 

SUB  GROUP 

18  SUBJECT  TERMS  (Continue  on  reverse  d  necessary  and  identity  by  block  number) 

ELF  Ionospheric  propagation 

VLI  Propagation  modeling 

Earth  ionosphere  waveguide 


1 9  ABSTRACT  ( Continue  on  reverse  d  necessary  and  identity  by  block  number) 

'  A  simple  surface  propagation  model  is  commonly  used  to  estimate  the  effects  of  ionospheric  disturbances  on  extremely  low 
frequency  (ELF)  propagation.  Though  often  adequate,  one  shortcoming  of  the  model  is  its  failure  to  allow  lor  excitation  by  the  broad¬ 
side  component  of  a  horizontal  dipole.  That  excitation  can  be  quite  significant  in  a  sporadic-E  environment,  w  hen  the  modal  polarization 
can  contain  a  substantial  mixture  of  TF.  component.  In  this  study  a  beginning  is  made  on  ELF  propagation  model  development  which 
allows  for  broadside  excitation  as  well  as  a  systematic  allowance  for  height  gain  effects.  The  development  is  for  a  flat  earth  with  a  laterally 
homogeneous,  cylindrically  symmetric  disturbance  centered  over  the  transmitter.  The  method  utilizes  normal-mode  decomposition. 
Matching  equations  at  the  boundary  of  the  disturbance  arc  developed  by  evaluating  height  gain  integrals,  and  it  is  shown  that  this  method 
is  tantamount  to  simply  matching  the  ground  vertical  electric  and  azimuthal  magnetic  fields  at  the  boundary.  Results  for  a  spoiadic-l 
type  environment  indicate  substantial  departures  from  the  surface  propagation  model  and  WKB  predictions. 
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I.  INTRODUCTION 


Since  an  ELF  signal  from  a  remote  transmitter  is  received  over  a  range  of 
azimuth  angles,  lateral  ionospheric  gradients  produced  by  sporadic-E  layering 
or  nuclear  depressions  can  produce  significant  effects  on  propagation  in  the 
lower  ELF  band.  This  is  because  the  Fresnel  zone  size  can  be  comparable  to 
the  transverse  dimensions  over  which  the  disturbed  ionosphere  changes 
significantly.  Although  a  number  of  workers  have  addressed  the  question  of 
propagation  in  waveguide  environments  which  vary  both  along  and  transverse  to 
the  path  of  propagation  [Wait,  1964;  Greifinger  and  Greifinger,  1977;  Field 
and  Joiner,  1979,  1982;  Pappert,  1980,  1985;  Ferguson,  Hitney  and  Pappert, 
1982;  Shellman,  1983],  no  formulation  exists  which  can  fully  account  for  the 
propagation  effects  produced  by  a  localized  disturbance  with  simultaneous 
allowance  for  vertical  inhomogeneity,  and  anisotropy  in  a  spherical 
geometry.  Thus,  it  has  been  common  practice  to  estimate  the  effects  produced 
by  localized  ionospheric  disturbances  by  using  a  simple  surface  propagation 
model  introduced  by  Wait  and  more  fully  developed  by  the  Grei fingers  and 
Field.  The  theory  is  predicated  on  the  assumption  that  the  field  can  be 
separated  into  lateral  and  height-dependent  factors,  and  the  model  has  been 
used  to  estimate  the  behavior  of  the  lateral  factor.  Excitation  factor  and 
height  gain  effects  are  generally  ignored.  When  applying  the  method  to 
nocturnal  environments,  additional  assumptions  are  made.  Among  these  is  the 
omission  of  the  TE  component  of  the  modal  polarization.  This  omission  is  well 
justified  in  the  ambient  case.  However,  it  is  known  that  under  sporadic-E 
layering,  considerable  mixing  between  TE  and  TM  waves  can  occur.  When  the 
surface  propagation  model  is  applied  to  nocturnal  sporadic-E  environments, 
this  TE  component  is  neglected. 


In  this  study  a  start  is  made  on  ELF  propagation  model  development  which 
allows  for  the  inclusion  of  excitation-factor  and  height  gain  effects  along 
with  allowance  for  the  TE  component.  The  development  is  for  a  flat  earth  with 
a  homogeneous,  cylindrically  symmetric  disturbance  centered  over  the 
transmitter  (see  Figure  1).  The  method  can  be  extended  to  more  complicated 
geometries  by  methods  discussed  by  the  Grei fingers  and  recently  implemented  by 
Pappert  [1985]  for  surface  model  calculations.  The  development  utilizes 
normal -mode  decomposition  (only  the  single  non-evanescent  mode  which 
propagates  at  ELF  has  been  used  for  the  calculations),  and  the  normal  mode 
parameters  are  taken  to  be  independent  of  the  direction  or  azimuth  of 
propagation.  A  method  akin  to  mode  conversion  has  been  used  to  develop  the 
matching  equations  at  the  boundary,  r0.  of  the  disturbance.  This  method  is 
shown  to  be  tantamount  to  simply  matching  the  ground  vertical  electric  and 
azimuthal  magnetic  fields  at  rQ.  Results  indicate  substantial  departures  from 
the  surface  propagation  model  and  WKB  predictions. 

Essential  theoretical  background  is  given  in  Section  II.  Azimuthal 
dependence  as  well  as  other  features  of  the  modes  are  discussed  in  Section  III 
along  with  the  final  field  results.  Conclusions  and  recommendations  are 
summarized  in  Section  IV. 
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II.  MODEL  AND  THEORY 


The  problem  is  illustrated  in  Figure  1,  where  a  horizontal  electric 
dipole  (current  moment  Idl)  is  shown  within  the  earth- ionosphere  waveguide  and 
beneath  the  center  of  a  large  cyl indrically  symmetric  patch  of  sporadic-E. 
The  vertical  composition  of  the  nighttime  ambient  electron  and  ion  profiles 
with  and  without  the  sporadic-E  layer  is  shown  in  Figure  2.  The  ground  is 
taken  to  be  homogeneous,  with  conductivity  Og  and  permittivity  eg.  A  crucial 
assumption,  discussed  more  fully  in  Section  III,  is  that  the  propagation  is 
independent  of  its  azimuth.  The  electric,  ft,  and  magnetic,  ft,  fields  excited 
by  the  horizontal  dipole  can  be  derived  from  two  scalar  potentials  V  and  W 
[Wait,  1963],  where  V  is  identified  with  quasi  TM  excitation  and  W  with  quasi 
TE  excitation.  Total  field  components  within  the  guide  are  given  by  [Webber 
and  Peden,  1971]. 
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where  w,  nQ,  cq  and  k  are  angular  frequency,  free  space  permeability,  free 
space  permittivity  and  free  space  wavenumber,  respectively.  A  time  dependence 
exp(j'dt)  is  assumed  and  j  = 

Free  space  field  components,  e£  ,  and,  HP,  generated  by  the  horizontal 
electric  dipole  are: 
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(3) 


rfiP  =  ~  Si n 4>  /  S2H,(2) 
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(kSr)exp(-jkC|z-z0l)de 


where 


M  =  Idl/(4iyu)Eo)  ,  S  =  si  no  ,  C  *  cosO  ,  n  =  (4) 

Gamma,  r,  is  the  contour  with  lower  endpoint  -n/2-j®  and  upper  endpoint 
n/2+j«.  h{2)  is  the  Hankel  function  of  order  1  of  the  second  kind.  The  plus 
sign  in  Eq.  (2)  applies  when  z  >  zQ  and  the  minus  sign  when  z  <  z0,  where  z0 
is  the  transmitter  height  (see  Figure  1). 

To  the  primary  field  components  represented  by  Eqs.  (2)  and  (3)  must  be 
added  the  field  components  associated  with  the  boundary  reflections  or, 
equivalently,  the  image  fields.  Let  R  be  the  plane  wave  reflection  matrix, 
referenced  to  ground  level,  associated  with  reflections  from  structure  above 
z0  and  R  the  corresponding  reflection  matrix  associated  with  reflections  from 
structure  below  zQ.  In  terms  of  the  elements  defined  by  Budden  [1961],  the 
reflection  matrices  associated  with  the  z  field  components  are: 


The  notation  i  denotes  vertical  polarization,  and  the  notation  l  denotes 
horizontal  polarization.  The  first  subscript  refers  to  the  polarization  of 
the  electric  field  of  the  upgoing  wave  while  the  second  subscript  refers  to 
the  corresponding  polarization  of  the  downcoming  wave.  Following  Pappert 
[1968],  the  total  z  field  components  within  the  guide  may  be  written  as 
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where 


F  =  (1-RR)'1 


(7) 


By  carrying  out  the  matrix  operations  indicated  in  Eqs.  (6)  and  (7)  and 
evaluating  the  integral  by  means  of  residue  theory  [Budden,  1962]  one  obtains 
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Here  m  is  a  mode  index  (only  one  non-evanescent  mode  exists  in  the  ELF  band) 
and  the  modal  function  is 


“  “-.R,  lRi>  -  x\  „R.  iR|  .R1 


(10) 


The  modal  equation  A  =  0  has  also  been  used  in  arriving  at  Eqs.  (8)  and  (9). 


From  Eqs.  (1-a)  the  total  potentials  V  and  W  within  the  guide  are 
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Equations  (11)  and  (12)  may  be  conveniently  written  as 
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and  the  modal  height  gains  hy  and  ey  are  defined  by 
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All  other  height  gains  are  given  by  Maxwell's  equations  subject  to  y 
invariance  and  an  exp(-jkSx)  dependence.  The  other  height  gains  are  (note 
that  h  is  defined  as  n  times  the  magnetic  intensity). 


6 


dh 


(17) 


e  -  -  i-  -X-  =  -C 
x  jk  dz 


e^Cz-  „R .  e 


-jkCz 


ae. 


1  +  uRn 


.  ejkCz  s  -jkCz 

,  h  -  i-  -X-  -  fC  - ii - 

jk  Sz  1  +  ±R ± 


e.  = 


jkexp(-jkSx)  dx 


-^-[exp(-jkSx)]  =  -Sh  ,  h  =  - 


y’  z 


jkexp(-jkSx)  5x 


—Cexp(-jkSx)]  = 


Sey  (18) 


Because  of  the  assumed  propagation  isotropy  in  azimuth,  the  height  gains  are 
independent  of  the  azimuth  angle  4>.  In  Eqs.  (16)  and  (17),  f  is  the  modal 
polarization  ratio  given  by 

f  =  ,»,)  =  ft+JJ  fiL  U9) 

hy(z-O)  ^1+I|R||^  1R II  1R1  (1+I|RI1^  1_1R1  iRl 

Thus  far  the  discussion  has  been  restricted  to  the  region  within  the 
guide.  The  modal  height  gains,  however,  have  meaning  both  within  and  wihout 
the  guide  (e.g.,  within  the  ionosphere)  and  are  calculable  by  means  of 
numerical  integration  of  Maxwell's  equations  [Pitteway,  1965;  Smith,  1970]. 
Therefore,  with  the  understanding  that  hy,  e^  ,  etc.,  represent  these  general 
height  gain  functions,  the  field  components  deduced  from  V  and  W  by  means  of 
Eqs.  (1)  can  be  continued  into  the  ionosphere. 

In  anticipation  of  solving  the  problem  posed  in  Fig.  1  by  satisfying 

continuity  requirements  at  the  boundary  r=rQ,  the  superscript  p  will  be  used 

to  denote  values  in  the  perturbed  region,  and  the  superscript  u  will  be  used 

to  denote  values  in  the  unperturbed  region.  To  the  primary  fields  in  the 

perturbed  region  must  be  added  the  secondary  fields.  The  crucial  assumption 
is  now  made  that  these  secondary  fields  generated  by  the  discontinuity  at 
r  =  r0  may  be  described  in  terms  of  a  normal-mode  expansion  so  that  the 
secondary  potentials  take  the  form 
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where  is  the  Bessel  function  of  the  first  kind  of  order  one  and  guarantees 
that  the  secondary  fields  are  finite  at  the  origin.  It  is  likewise  assumed 
that  the  potentials  in  the  unperturbed  region  can  be  written  in  terms  of  a 
normal -mode  expansion  so  that 
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IV  SUMMARY 


In  this  report,  a  development  to  allow  for  the  inclusion  of  excitation 
factor  and  height  gain  effects  on  ELF  propagation  in  a  laterally  inhomogeneous 
guide  has  been  started.  The  development  is  for  a  flat  earth  geometry  and  the 
simple  environment  of  Figure  1.  Several  crucial  assumptions  have  been  made. 
First,  it  has  been  assumed  that  the  perturbing  environment  is  homogeneous  as 
well  as  cyl indrically  symmetric.  Second,  it  has  been  assumed  that  the  primary 
fields,  as  well  as  secondary  fields  generated  by  the  boundary  at  r  =  r  ,  are 
described  by  normal  modes,  and  only  the  consequences  of  the  single  non- 
evanescent  mode  have  been  investigated.  Third,  the  normal- mode  parameters 
have  been  taken  to  be  independent  of  the  direction  or  azimuth  of  propagation. 

A  method  akin  to  mode  conversion,  termed  in  this  study  the  SMA  (for 
Single  Mode  Analysis)  method,  has  been  used  for  establishing  the  boundary 
equations  at  r  =  rQ.  As  in  mode  conversion  analysis,  the  SMA  method  involves 
numerical  integration  of  height  gain  functions.  It  has  been  shown  that  the 
SMA  method  is  tantamount  to  matching  H  and  Ez  at  r  =  rQ  at  the  ground.  This 
is  a  significant  result  because  it  implies  the  reasonableness  of  a  generalized 
surface  propagation  model  which  allows  in  a  systematic  way  for  excitation 
factor  and  height  gain  effects  as  well  as  broadside  excitation. 

It  has  been  found  that  excitation  factor  and  height  gain  effects  can  be 
quite  significant  for  sporadic-E  environments  conforming  to  the  geometry  of 
Figure  1.  In  particular,  broadside  excitation,  because  of  the  sizable  TE 
contamination  of  the  modal  polarization,  can  be  very  substantial.  TE 
contamination  would  also  be  expected  to  play  a  role  were  the  disturbance  to 
fall  over  the  receiver  rather  than  over  the  transmitter. 
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the  dominating  terms  of  Eqs.  (29)  through  (32)  involve  integrals  of  products 
of  ez  and  h^  terms.  These  terms  are  nearly  constant  and  of  order  unity  below 
about  73  km  (see  Figures  29  and  30),  where  most  of  the  contribution  to  the 
integrals  comes  from.  The  integrals  involving  the  products  of  TE  components 
would  generally  be  much  smaller,  as  Figures  29  and  30  show.  At  any  rate,  it 

is  a  fact  that  simply  matching  the  and  Ez  components  at  r  =  r0  at  the 

ground  is  for  all  practical  purposes  equivalent  to  the  height  gain  integration 
method  (the  SMA  method)  used  in  Section  U  to  develop  the  results.  This  is 
made  clear  in  Figures  31  through  34,  which  show  a  worst-case  comparison 
between  SMA  results  and  "Field  Matching"  results  (i.e.,  results  obtained  by 

demanding  continuity  of  and  Ez  at  the  ground  at  r  =  rQ)-  The  differences 
are  indeed  very  slight.  This  is  a  significant  result  because  it  implies  the 

reasonaDl eness  of  a  generalized  surface  propagation  model  which  allows  in  a 
systematic  way  for  excitation  factor  and  height  gain  effects  as  well  as 
broadside  excitation. 

A  disquieting  though  not  surprising  feature  of  the  analysis  is  that  the 
TE  components  are  not  continuous  at  r  =  rQ.  This  would  be  expected  for 
the  e  ^  and  Hz  components  since,  as  Figures  29  and  30  show,  ey  and  hz  vary 

substantially  across  the  guide.  However,  the  hx  component  is  quite  sub¬ 

stantial  and  reasonab1y  stable  across  the  guide.  But  as  Figure  35  (a  repre¬ 
sentative  example)  shows,  the  Hr  component  experiences  a  large  discontinuity 
at  r  -  rQ  at  the  ground.  This  is  not  surprising  because,  as  we  have  just 
seen,  the  SMA  method  is  equivalent  to  matching  only  the  H ^  and  Ez  components 
at  r  =  r-j  at  the  ground.  It  is  not  clear  what  the  resolution  to  this  dilemma 
is.  Perhaps  inclusion  of  evanescent  modes  would  improve  the  agreement  for 

continuity  of  Hr  without  destroying  the  and  Ez  continuity.  However,  no 

efforts  along  those  lines  have  been  attempted  in  the  current  study. 
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close  to  broadside,  excitation  variation  for  the  unperturbed  guide  is  close  to 

the  surface  model  result;  however,  the  perturbed  guide's  normalized  source 

height  gain  variation  shows  marked  departure  from  the  surface  model  result. 

This  is  because  of  the  large  TE  admixture  to  the  modal  polarization  in  the 

perturbed  region.  Whereas  for  the  unperturbed  guide  I  eu(o) I >> I eu ( o ) I ,  the 

x  y 

situation  for  the  perturbed  guide  is  that  | e P ( o ) I  is  somewhat  greater 
than  |eP(o)|,  and  that  results  in  the  interference  null  in  Figure  16  occurring 
at  —  36°  and  for  the  region  $  <0  being  essentially  a  region  of  constructive 
interference. 

Figures  17  and  18  show  results  for  the  amplitude  radios  IhL/H^I  and 

4>  4> 

IEZ/E^|  at  the  ground  for  crossed  dipole  sources  which  are  out  of  phase  by 
90°.  Shown  are  SMA,  WKB  and  surface  propagation  model  results.  Again  the  SMA 
results  are  very  nearly  continuous  at  r  =  rQ  while  the  WKB  results  manifest 
sizable  discontinuities  there.  Exterior  to  the  disturbance,  the  SMA  and 
surface  propagation  model  results  are  in  reasonable  agreement.  Interior  to 
the  disturbance,  the  SMA  results  for  differ  significantly  from  the  surface 
propagation  model  results  while  the  interior  results  for  E2  for  the  two  models 
are  in  reasonable  agreement. 

Figures  19  through  28  show  SMA  and  WKB  phase  results  for  and  Ez  for 
the  cases  for  which  amplitude  results  have  been  given  previously.  The  signi¬ 
ficant  features  are  the  near  continuity  of  the  phase  of  h  and  E,  at  r  =  rn 

<p  £  o 

for  the  SMA  method  and  the  large  discontinuity  exhibited  there  by  the  WKB 
results.  The  near  continuity  of  the  amplitudes  and  phases  at  r  =  rQ  means 
that  the  height  gain  integration  method  used  in  Section  II  to  develop  the 
results  is  tantamount  to  simply  matching  H  and  Ez  components  at  r  =  r0  at  the 
ground.  The  latter  method  involves  no  height  gain  integrations  such  as  used 
in  the  method  of  Section  II.  It  is  suspected  that  the  reason  for  this  is  that 


iii)  SMA  results  show  significant  departures  from  both  the  surface 
propagation  model  and  the  WKB  results. 

iv)  At  short  ranges  the  SMA  and  l*B  results  generally  show  good 
agreement. 

v)  In  contrast  to  the  surface  propagation  model  results,  the  SMA  results 
show  a  marked  angular  (i.e. ,  4>)  dependence. 

Result  (i)  above  is  due  in  large  measure  to  the  large  discrepancy  between 
the  perturbed  ( 9p=  57.16°  -  65.03°j)  and  unperturbed  (9^  =  84.48  -  34.12  j) 

eigenangles  as  well  as  the  difference  between  the  perturbed  [ep(o)  =  (5.78 

-  6.77j)xl0-4]  and  unperturbed  [ej(o)  =  (5.25-2.21j)xl(T5]  field  components 
and  the  difference  between  the  excitation  factors  ( \p  s  2.32  -  4.99j  and  \u  = 

-  4.32xl0"2  -  7 .42 j ) .  In  connection  with  (ii)  above,  it  is  pointed  out  that 
continuity  of  the  field  components  ar,d  Ez  is  not  obviously  forced  by  the 
matching  method  of  Section  II.  Departures  mentioned  in  (iii)  point  out 
clearly  the  significance  of  height  gain  and  excitation  factor  effects  ignored 
in  the  surface  propagation  model.  Observation  (iv)  is  explained  by  the  fact 

that  the  secondary  fields  [which  depend  upon  J^(kSPr)  and  Jj(kSPr)]  for 

I kSPr I <<1  are  small  compared  with  the  primary  fields  which  depend 
upon  n[2)(kSpr)  and  H[2)’(kSpr).  Figure  16  illustrates  the  source  of  the 

angular  dependence  mentioned  in  (v).  Shown  is  the  normalized  source  height 
gain  (e  (o)cos*  +  e  (o)sino)/e  (o)  in  dB  for  both  the  perturbed  "p"  and 

unperturbed  "u"  quantities  versus  azimuth.  Shown  also  is  the  201  og^Q  (cos4>( 
variation  predicted  by  the  surface  propagation  model  (which  assumes  only 
excitation  from  the  end-on  component  of  the  source  dipole).  Except  for  angles 


vertical  component  of  the  geomagnetic  field  for  high  latitudes.  Azimuthal 
invariance  is  strictly  valid  for  a  dip  angle  of  90°. 

All  remaining  figures  in  this  study  are  for  75  Hz,  and  for  a  vertical 

geomagnetic  field  (i.e.,  dip  of  90°)  of  strength  .41  x  sin(77°)  gauss.  The 

ground  conductivity  has  been  taken  to  be  10"2  siemens/m  and  the  dielectric 

constant  15.  The  lower  conductivity  of  about  2.8  x  10“^  siemens/m  at  the 

transmitter  can  be  included  by  simply  multiplying  the  height 

gains  ej(p)  and  e£(o)  by\10“2/2.8  x  10“4  =  5.98  (or  15.5  dB).  However,  the 
x  y 

factor  is  of  no  consequence  in  the  present  study,  where  either  amplitude 
ratios  or  phase  differences  are  plotted.  Finally,  the  radius  of  the 
disturbance  has  been  taken  to  be  500  km. 

Figures  8  through  15  show  results  for  the  amplitude  ratios  |H^/H^|  and 
|EZ/E“|  at  the  ground  for  an  electric  dipole  source  at  the  ground  oriented  in 
the  x  direction  (see  Figure  1).  Figures  8  and  9  are  for  broadside  launch 
(i.e.,  <t>  =  +90°),  Figures  10  and  11  are  for  $  =  45°,  Figures  12  and  13  are  for 
end-on  launch  (i.e.,  <j>  =  0°),  and  Figures  14  and  15  are  for  <t>  =  -45°.  Results 
of  the  present  study  are  labeled  SMA  for  Single  Mode  Analysis.  Shown  also  on 
Figures  10  through  15  are  results  and  results  of  the  surface  propagation 
model,  which  neglects  excitation  factor  and  height  gain  effects  (Pappert, 
1985).  The  surface  propagation  results  are  $  independent.  There  are  several 
observations  that  can  be  made  from  these  figures. 

i)  WKB  results  generally  show  a  sizable  discontinuity  at  r  =  rQ. 

ii)  SMA  results  are  very  nearly  continuous  at  r  =  rQ  though  range 
derivatives  are  discontinuous. 
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ELF  propagation  in  the  ambient  guide  can  be  taken  to  be  independent  of 
azimuth. 

Figure  5  shows  that  the  azimuthal  deviation  from  the  average  attenuation 
rate  for  the  disturbed  profile  (i.e.,  ambient  plus  sporadic-E)  is  less  than 
10%.  The  large  difference  between  the  mean  attenuation  rates  for  the  two  dips 
shown  on  Figure  5  results  from  the  critical  role  which  the  vertical  component 
of  the  geomagnetic  field  plays  in  determining  the  attenuation  rate.  This  is 
discussed  in  the  work  of  Pappert  and  Moler  [1978].  Possible  contributions  of 
this  effect  to  azimuthal  variations  as  well  as  azimuthal  dependence  due  to 
possiDle  lateral  inhomogneity  of  the  layer  are  ignored  in  this  study. 
Finally,  Figure  6  shows  that  the  azimuthal  deviation  from  the  average 
normalized  phase  velocity  for  the  disturbed  profile  is  also  less  than  10%. 
Thus,  though  the  evidence  for  assuming  azimuthally  independent  propagation  in 

this  instance  is  not  as  formidable  as  in  the  case  of  the  ambient  profile,  the 
azimuthal  variations  for  constant  dip  angle  are  still  sufficiently  small  to 
make  reasonaDle,  at  least  as  a  first  approximation,  the  assumption  of 
azimuthal  invariance.  The  assumption  is,  of  course,  most  reasonable  when 
applied  to  effects  caused  by  natural  or  artificial  ionospheric  depressions  and 
has  served  as  the  basis  for  the  theoretical  development  of  Section  II. 

Figure  7  shows  attenuation  rate  as  a  function  of  frequency  for  the 
ambient  nighttime  profile  and  the  ambient  plus  sporadic-E  profile.  Two  curves 
for  the  disturbed  environment  are  shown.  The  solid  curve  is  for  a  dip  angle 

of  77°,  a  geomagnetic  field  strength  of  0.41  gauss,  a  ground  conductivity  of 
1Q~2  siemens/m  and  dielectric  constant  of  15.  The  dashed  curve  is  for  a  dip 

of  90°  and  a  geomagnetic  field  strength  of  0.41  x  sin(77°)  gauss  (i.e.,  the 
radial  component  of  the  geomagnetic  field  used  for  the  solid  curve).  The  near 
coincidence  of  the  two  curves  demonstrates  the  controlling  feature  of  the 
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III.  RESULTS 


In  a  previous  study  [Pappert  and  Moler,  1978]  full-wave  outputs  for  an 
ambient  nighttime  profile  disturbed  by  a  sporadic-E  layer  (see  Figure  2)  have 
been  analyzed  with  respect  to  ionospheric  absorption  and  reflection  features 
at  lower  ELF.  In  that  work  it  was  shown  that  an  order- of-magnitude 
enhancement  1r.  the  attenuation  rate  due  to  the  sporadic-E  layering  was  a 
distinct  possibility  (see  also  Figures  3,  5  and  7  herein).  In  this  section, 
features  of  the  propagation  associated  with  the  simple  sporadic-E  environment 
of  Figures  1  and  2  will  be  explored  by  implementation  of  theory  developed  in 
Section  II. 

Figures  3  through  6  show  the  behavior  of  the  attenuation  rate  and  phase 
velocity  (normalized  to  the  free  space  speed  of  light)  at  75  Hz  as  a  function 
of  azimuth  of  propagation  relative  to  magnetic  north  for  dip  angles  of  60°  and 
75°.  A  conductivity  of  a  =  10'^  siemens/m,  dielectric  constant  of  15,  and 
geomagnetic  field  strength  of  0.41  gauss  have  been  used  for  the 
calculations.  Figure  3  indicates  mid  and  northerly  latitude  ambient 

attenuation  rates  of  slightly  greater  than  0.8  dB/Mm.  This  is  somewhat  lower 
than  the  WTF  measurements,  which  estimate  a  nighttime  mid-latitude  value  of 
1.0  dB/Mm.  This  comparison  could  be  indicative  of  a  deficiency  in  the 
ionospheric  model,  although  other  models  also  generally  yield  lower 
attenuation  rates  than  those  suggested  by  WTF  measurements.  Most  measurements 
of  ELF  nocturnal  fades  have  been  made  for  paths  for  which  the  dip  angle 
is  >  60°.  The  curves  of  Figure  3  indicate  an  azimuthal  variation  in 
attenuation  rate  of  less  than  0.03  dB/Mm.  Similarly,  Figure  4  shows  the 
difference  In  the  normalized  phase  velocity  for  the  east  to  west  and  west  to 
east  directions  is  less  than  1%.  Thus,  for  paths  of  concern  in  this  study. 
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As  a  matter  of  interest,  comparisons  will  also  be  made  in  the  following 
section  with  the  so-called  WKB  result  [Bickel  et  al.,  1970]  which  involves  the 
geometric  mean  of  the  excitation  factor  at  the  terminals  of  the  path  and  the 
average  of  the  propagation  constant  over  the  path.  The  WKB  formulas  for  the 
field  as  implemented  in  the  present  study  are  as  follows: 
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Since  ratios  of  the  perturbed  to  unperturbed  fields  are  the  quantities 
displayed  in  the  following  section,  for  convenience  the  unperturbed  fields  are 
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Equations  (29)  through  (32)  determine  (xj,  x2,  x3,  x4)m  and  in  terms  of 
these  quantities  the  field  components  discussed  in  the  following  section  are 
developed  from  Eqs.  (1),  (13),  (14)  and  (20)  through  (23).  They  are 
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In  Eq.  (28)  the  dagger  sign,  t,  indicates  the  adjoint  (i.e.,  the  complex 
conjugate  transpose). 

To  generate  the  field  matching  equations,  first  multiply  Eq.  (24)  through 
by  cos <j>  and  integrate  $  from  0  to  2 it.  Next  multiply  through  by  [qP(- ) and 
integrate  over  z.  If  N  is  the  total  number  of  modes  used,  this  generates  one 
equation  for  the  4  x  N  unknowns  (x^,  xg,  X3,  X4 ) m,  where  m  takes  on  the  values 
1  through  N.  Next  multiply  the  equation  obtained  from  the  $  integration 
by  (qp(+) )  and  integrate  over  z.  This  generates  a  second  equation.  By 
extracting  the  sin$  terms,  two  more  equations  are  generated  in  a  similar 
fashion.  Performing  these  operations  for  n  =  1  through  N  generates  the 
requisite  4  x  N  equations.  The  equations  are  (here  the  mode  indices  are 
included  explicitly): 
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As  discussed  by  Pappert  and  Smith  [1972],  there  exists  height  gains 
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associated  with  an  adjoint  guide  for  which  the  inner  products  [qp(-),  fp(-)] 
and  [qJJ(+),  f^(+)]  are  zero  when  m  ^  n.  The  matching  equations  will  be 
developed  by  a  method  somewhat  akin  to  the  treatment  of  mode  conversion  in  the 
VLF  band  [Pappert  and  Snyder,  1972],  although  it  is  by  no  means  clear  in  the 
present  case  that  it  is  an  optimum  or  preferred  procedure.  As  will  be  seen 
later,  the  method  does  have  the  convenient  consequence  that  it  is  essentially 
equivalent  to  simply  matching  and  Ez  at  r  *  r0  at  the  ground.  As  a 
preliminary  to  developing  the  matching  equations,  the  following  definitions 
will  be  useful: 


jH[2)'(kSpro)eP(z) 

jJ1(kSprQ)ep(z) 

eP(z) 

pm  “ 

H}2)(kSPro)eP{z) 

jH}2)'(kSpro)hp(z) 

»  ^m  = 

J i (kSpr^)ep(  z) 
jjj(kSpr0)hp(z) 

*  Rm 

0 

t>P(2) 

h|2) (kSprQ)hp(z) 

m 

J,  (kSpr  )hp(z) 

^  1  U  L 

m 

0 

jh{2)  (kSurQ)ey{z) 

”ej(zf 

nj2^(kSUr0)e^(z) 

0 

^m  = 

jHj2)’(kSur0)hJ(z) 

.  um  = 

hUx<z) 

H{2)(kSurQ)h“(z) 

m 

_  0 

J,(kSpr  Q) 

+jSp - - -  ( -x,  si n<t>  +  x9cos<*>) 

kS pr  1  £ 


-  -k2 X<!su ( cos <t>  +  x4 si n<j>) 
m 


eP(2) 


hP(2) 


jH[2)’(kSur0)eJ(2) 

H[2)(kSur0)ep(2) 

jH[2)'(kSur0)hp(2) 

h|2) (kSur  )hp(z) 
l  02 


(Z)  u  ex(z) 

i.  H i  (kSur  )  * 

+  jSu - - -  ( -XoSi n<t>  +  x.cos$)  (24) 

kSur„  3  hJU) 

0 

*-  -•  m 

In  Eq.  (24)  the  top  element  of  the  column  vector  is  E^  ,  the  second  Ez»  the 

third  rrl  and  the  last  is  The  left-hand  side  of  the  equation  gives  the 

w  z 

fields  associated  with  the  perturbed  region  at  r  =  rQ  and  the  right-hand  side 
gives  the  fields  associated  with  the  unperturbed  region  at  r  =  r0* 

With  a  finite  number  of  modes  (and  in  the  present  case  only  one  non- 
evanescent  mode  exists)  it  is  impossible  to  satisfy  Eq.  (24)  for  all  z.  Thus, 
the  continuity  will  be  developed  only  approximately  following  a  prescription 
suggested  by  asymptotic  considerations  when  | kSPr0 I >>1  and  | kSur0 I >>1 .  In 
that  instance  the  leading  terms  of  the  forward  (+)  and  backward  (-) 
propagating  waves  associated  with  Eq.  (24)  becomes  a  product  of  radial 
functions  and  the  column  vectors 


*5<*>  ■' 


*ej(2) 

ej(z) 

±hj(z) 

hP(z) 


• 


e$U) 

ej(z) 

h“(z) 


A  disquieting  feature  of  the  analysis  is  that  the  TE  components,  and  in 


particular  the  Hr  component,  are  not  continuous  at  r  =  rQ.  This  is  not 
terribly  surprising  because,  as  just  mentioned,  the  method  used  to  obtain  the 

boundary  equations  is  equivalent  to  matching  only  the  and  Ez  components  at 
r  =  r0  at  the  ground.  Perhaps  inclusion  of  evanescent  modes  would  improve 
continuity  of  Hr  without  destroying  the  and  Ez  continuity.  However,  no 
effort  along  those  lines  has  been  attempted  in  the  current  study. 

Several  possible  extensions  of  the  present  effort  are: 

i)  Examine  effects  of  non-evanescent  modes  on  continuity  of  major  field 
components. 

ii)  Extend  analysis  to  allow  for  general  location  of  the  disturbance 
relative  to  the  transmitter. 

iii)  Extend  analysis  to  allow  for  radial  and/or  azimuthal  variations. 

iv)  Include  the  effect  of  earth  curvature. 

v)  Extend  the  analysis  to  el  1 i pti cal ly  symmetric  disturbances. 
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gure  1.  Source  and  waveguide  geometry. 
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Figure  22.  Disturbed  and  undisturbed  phase  difference  for 
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Figure  21.  Disturbed  and  undisturbed  phase  difference 
for  Ha  vs.  range. 
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Figure  20.  Disturbed  and  undisturbed  phase  difference  for 
Ez  vs.  range. 
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Figure  19.  Disturbed  and  undisturbed  phase  difference 
for  H  .  vs.  range. 
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Figure  17.  Ratio  of  disturbed  to  undisturDed  H  in  dB  vs. 
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Figure  16.  Excitation  efficiency  of  height  gains  vs 
azimuth. 
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Figure  13.  Ratio  of  disturbed  to  undisturbed  Ez  in  dB  vs. 
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Disturbed  and  undisturbed  phase  difference 
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Figure  32.  Comparison  of  SMA  and  Field  Matching  phases 
for  iz. 


60 


0  100 
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Figure  34.  Comparison  of  SMA  and  Field  Matching  phases 
for  Ez. 
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Figure  35.  Ratio  of  disturbed  to  undisturbed  Hp  in  dB  vs. 
range. 
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